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ABSTRACT 

We present a three-dimensional density model of coronal prominence cavities, 
and a morphological fit that has been tightly constrained by a uniquely well- 
observed cavity. Observations were obtained as part of an International Helio- 
physical Year campaign by instruments from a variety of space- and ground-based 
observatories, spanning wavelengths from radio to soft-X-ray to integrated white 
light. From these data it is clear that the prominence cavity is the limb man- 
ifestation of a longitudinally-extended polar-crown filament channel, and that 
the cavity is a region of low density relative to the surrounding corona. As a 
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first step towards quantifying density and temperature from campaign spectro- 
scopic data, we establish the three-dimensional morphology of the cavity. This is 
critical for taking line-of-sight projection effects into account, since cavities are 
not localized in the plane of the sky and the corona is optically thin. We have 
augmented a global coronal streamer model to include a tunnel-like cavity with 
elliptical cross-section and a Gaussian variation of height along the tunnel length. 
We have developed a semi- automated routine that fits ellipses to cross-sections 
of the cavity as it rotates past the solar limb, and have applied it to Extreme 
Ultraviolet Imager (EUVI) observations from the two Solar Terrestrial Relations 
Observatory (STEREO) spacecraft. This defines the morphological parameters 
of our model, from which we reproduce forward-modeled cavity observables. We 
find that cavity morphology and orientation, in combination with the viewpoints 
of the observing spacecraft, explains the observed variation in cavity visibility for 
the east vs. west limbs. 

Subject headings: solar corona, magnetic flux ropes, coronal mass ejections 


1. Introduction 

Coronal mass ejections (CMEs) are spectacular solar eruptions and the primary drivers 
of “space weather” at the Earth (Pick et al. 2006; Hudson et al. 2006). CMEs are most 
commonly observed in white light, and many possess the classic “three-part” morphology 
of a bright expanding loop, followed by a relatively dark cavity, and lastly a bright core 
associated with an erupting prominence. Such three-part CMEs have been observed to 
erupt from previously quiescent cavities with embedded prominences (Sterling & Moore 
2004; Gibson et al. 2006). Indeed, prominences and cavities are ubiquitous features of the 
non-erupting corona: several may be visible on any given day. In the polar crown they can be 
long-lived (on the order of months), either reforming after or only partly erupting in CMEs 
(Gibson &; Fan 2006; Su et al. 2006; Liu et al. 2007; Tripathi et al. 2007, 2008). Prominences 
and their cavities thus provide clues to magnetohydrodynamic (MHD) equilibrium states of 
the solar corona, and to the processes which can destabilize these equilibria and drive certain 
coronal mass ejections. The thermal and dynamic properties of cavities provide information 
about their magnetic structure and MHD equilibrium (Zhang & Low 2004; van Ballegooijen 
&; Cranmer 2008). It is therefore critical to obtain information about cavity plasma density, 
velocity, and temperature via multi-wavelength observations in order to test and constrain 
MHD models. 

Cavity densities are most unambiguously determined from white-light observations, 
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which are temperature- independent, and these indicate that cavity density has a lower limit 
of approximately half the density of a surrounding streamer at the same height (Fuller et al. 
2008; Fuller &; Gibson 2009). Radio observations have been consistent with this (Marque 
(2004) and references therein). 

Studies of velocities associated with the non-erupting cavity have more recently been 
undertaken: Schmit et al. (2009) observed large-scale flows on the order of 5-10 km/sec 
in coronal cavities using extreme ultraviolet (EUV) and infrared spectroscopy. Spinning 
motions within cavities have also been observed, with a preferred direction of spin consistent 
with meridional flows (near solar minimum) (Wang & Stenborg 2010). 

Temperature determination of cavities is a more challenging problem. From white light, 
the only information that can be obtained is an inferred density scale height and associated 
“hydrostatic temperature” (Guhathakurta et al. 1992). A comparison of white- light scale 
heights in cavities vs. their surrounding streamers indicates a higher hydrostatic temperature 
for most cavities (Fuller et al. 2008; Fuller & Gibson 2009). Since hydrostatic balance is along 
magnetic flux tubes, however, this implicitly assumes a uniform lower boundary density 
for all flux tubes that intersect the radial cuts along which scale heights are determined. 
This assumption may be particularly unsuitable for prominences and their cavities: indeed, 
an isothermal, MHD simulation of an emerging magnetic flux rope showed that density 
variations at the lower boundary led to differing cavity vs. streamer density scale heights as 
observed in white light, with no need for enhancement of cavity temperature (Fuller et al. 
2008). 

Independent information about cavity temperature can be obtained from observations 
of emission due to coronal spectral lines. Guhathakurta et al. (1992) examined densities and 
temperatures in a coronal cavity using eclipse observations in white light and the coronal red 
(6374 A Fex) and green (5303 A Fe XIV) lines. The white-light data indicated a clear density 
depletion, and, as in the Fuller et al. (2008) analysis, a hydrostatic temperature hotter than 
the streamer. On the other hand, their analysis of the line ratio of the red and green lines, a 
measurement which does not make hydrostatic assumptions, indicated cooler temperatures 
in the cavity. More recently, Habbal et al. (2010) used eclipse observations in white light 
and a range of iron lines to argue that low-emission regions around prominences were hot 
relative to the surrounding corona. However, the white-light observations associated with 
these structures indicated little or no density depletion. Thus, as the authors pointed out, 
they were not good examples of “cavities” in the sense of a lack of density, despite the fact 
that there was low emission at some wavelengths surrounding the prominence. Because the 
corona is optically thin, the orientation and three-dimensional morphology of the cavity may 
result in a large fraction of the line of sight being filled with non-cavity plasma, impacting 
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both density-depletion and temperature diagnostics. Thus it is important to consider the 
three-dimensional morphology of any given cavity to understand the line-of-sight projection 
effects. 

One approach for dealing with such projection effects is tomographic reconstruction, 
which requires no assumptions about boundary conditions or cavity morphology and ex- 
plicitly accounts for line-of-sight effects. Vasquez et al. (2009, 2010) used STEREO EUVI 
images to tomographically reconstruct coronal cavities, and found Differential Emission Mea- 
sure (DEM) consistent with lower density and higher temperature in the cavity than in the 
surrounding corona. This is the strongest evidence for “hot” cavities to date (in cavities 
which are demonstrably low-density), although it must be noted that these analyses have so 
far been limited to three EUVI bandpasses which cannot define a complicated DEM (Schmelz 
et al. 2007). It is possible that cavity temperature may vary internally, with only sub-regions 
being particularly hot - for example, near the cavity center where soft X-ray emission has 
been observed surrounded by an otherwise low-emission cavity (Hudson et al. 1999; Hudson 
& Schwenn 2000). (The “hot shrouds” around prominences observed by Habbal et al. (2010) 
could also have their origin in such structures). Observations with broad spectral range, 
plane-of-sky spatial resolution, and information about the cavity in three dimensions are 
required to conclusively establish cavity temperature relative to the surrounding corona. 

In this paper, we present observations of a coronal cavity observed in August 2007 during 
a multi-instrument observing campaign organized under the auspices of the International 
Heliophysical Year (IHY). These observations represent the most comprehensive observations 
of a cavity yet taken, at wavelengths ranging from radio to soft X-ray and with cavity- 
specific spectroscopic density and temperature diagnostics. Our ultimate goal is to use these 
diagnostics and observations of the full range of wavelengths to constrain a forward model of 
the three-dimensional density and temperature of the cavity. As a first step, in this paper we 
will determine the cavity’s three-dimensional morphology (eight free parameters, described 
below), using observations taken from a range of wavelengths and viewing angles. In Section 
2 we will describe the data taken during the IHY campaign, and demonstrate that our cavity 
is a region of low density and that it extends in longitude. In Section 3 we present our model 
and fit it to the data. In Section 4 we use the results of the forward model to examine how 
integration along the line of sight may produce differences in cavity visibility for different 
viewpoints. In Section 5 we present our conclusions. 
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Fig. 1. — 195 A images of cavity from STEREO EUVI-A, (a-f) Aug. 7-12, 2007. Numbers below 
the time stamp indicate Carrington longitude at the east limb on that day (90° less than the central 
meridian Carrington longitude). 

2. IHY Campaign on Coronal Prominence Cavities, August 8-18, 2007 

As described above, quantitative studies of density and temperature have to date pri- 
marily been done using radio and white-light observations. Cavities have also been observed 
in the emission corona at EUV and soft X-ray (SXR) wavelengths (Hudson et al. 1999; 
Hudson &; Schwenn 2000; Sterling & Moore 2004; Heinzel et al. 2008). To our knowledge, 
the IHY campaign was the first time a campaign was run with targetted observing plans 
customized for cavity density /temperature spectral diagnostics, with long integration times 
and a broad spectral range. It is essential to consider the three-dimensional morphology 
and orientation of cavities, to properly account for projection effects of non-cavity material 
into the line-of-sight (Wiik et al. 1994; Fuller et al. 2008). A critical campaign goal thus 
was to obtain multiple views and observations at a broad range of wavelengths, in order to 
constrain the three-dimensional morphology of the cavity and its surrounding streamer. 

Campaign observations included white-light polarized Brightness (pB) data from Mauna 
Loa Solar Observatory (MLSO) Mk4 K-coronameter, SXR data from the Geostationary Op- 
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erational Environmental Satellite (GOES) Solar X-Ray Imager (SXI), EUV image data from 
STEREO Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) EUVI A 
and B, Solar and Heliospheric Observatory (SOHO) Extreme ultraviolet Imaging Telescope 
(EIT) and Transition Region and Coronal Explorer (TRACE), UV and EUV spectroscopic 
data from SOHO Coronal Diagnostic Spectrometer (CDS) and Hinode Extreme ultraviolet 
Imaging Spectrometer (EIS), and radio observations from the Nangay Radioheliograph. We 
have also made use of H-a synoptic charts to determine the location of filament channels 
and magnetic neutral lines on the disk (McIntosh 2003). 

We observed portions of a well-established and longitudinally-extended polar crown 
filament channel (PCFC) over the period August 7-18, 2007. Targetted observations took 
place at the northeast limb (August 8-13), central meridian (August 14), and at the northwest 
limb (August 15-18). The portions of the PCFC we observed contained a prominence at least 
part of the time, and was clearly associated with a cavity at the northeast limb (see Figure 
1). For reasons we will discuss in Section 4, the cavity was not clearly visible when it 
rotated around to the northwest limb. Our analysis therefore focuses on fitting the model 
to observations of the cavity at the northeast limb. 
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Fig. 2. — Cavity on Aug. 9, 2007 at a range of wavelengths: (a) SOHO EIT 284 A , 01:04:50 UT 
(greyscale and white contours); (b) TRACE 171 A , 00:01:55 UT (inner image greyscale and red 
contours) embedded in EIT 284 A (outer image greyscale and white contours); (c) GOES-SXI, 
00:03:06 UT (greyscale and yellow contours); (d) EIT 284 A (greyscale) with all three contours 
overplotted. 


2.1. Cavity variation with wavelength 

Figure 2 shows observations of the cavity on August 9, 2007 near 00:00:00 UT, when it 
was clearly visible at the NE limb. Although there is some variation of bright substructures 
within the cavity, the basic size and shape of the cavity does not vary much with wavelength. 
Observations shown are from similar times, but range in wavelength and hence plasma- 
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Fig. 3.— (a) TRACE 171 (Aug. 9, 19:46:20 UT), and (b) SXI (Aug. 10, 00:03:06 UT) observations 
of the cavity, (c) TRACE 171 intensity vs. colatitude (north pole is to the right), and d) SXI 
intensity vs. colatitude. 

temperature sensitivity. The EUV 171 A observations have a peak temperature response 
at 1MK, and the 284 A observations have a peak temperature response at 2.2 MK (Dere 
et al. 2000; Howard et al. 2008). GOES-SXI has a broad temperature response at SXR 
wavelengths, and is sensitive to plasmas of temperatures greater than ~ 1 or 2 MK, with 
the sensitivity increasing somewhat with temperature up to ~ 10 MK (Pizzo et al. 2005). 

Figure 3 shows data from approximately 24 hours later. Three intensity vs. colatitude 
profiles are shown for a sequence of radial heights, illustrating that the cooler, 171 A emission 
and the hotter, SXR emission have similarly located cavity edges - essentially the inflection 
point of the intensity vs. colatitude profile (see Gibson et al. (2006) for discussion). This 
indicates that the intensity gradient at the cavity outer boundary is likely to be a density, 
rather than a temperature effect. Near the embedded prominence, however, the EUV lines 
must be interpreted carefully. First of all, a lack of emission could be due to either a 
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Fig. 4. — Cavity and associated prominence observed by EIT on August 10, 2007: (a) 284 A (01:06 
UT) and (b) 304 A (01:19 UT) 

density depletion or to absorption by dense, relatively cool prominence material. Secondly, 
transition-region and coronal emission may occur in the vicinity of a prominence; for example, 
Kucera &; Landi (2006) reported emission in spectral lines produced at LogT = 5.8 associated 
with an activated prominence (see also Labrosse et al. (2010)). This temperature is near 
the peak response of the 171 A bandpass, and so this may explain the central bump in 
intensity vs. colatitude in Figure 3 (a) and (c). Figure 4 (a) also shows 284 A emission 
colocated with the central prominence, and this may also in part be due to the presence 
of some transition-region lines in that bandpass (see Figure 12 in Dere et al. (2000) and 
Figure 7 in Howard et al. (2008)). Alternatively, the 284 A emission may be related to 
so-called “chewy nougats” , large regions of hot material (approaching « 2 MK) interior to 
cavities observed in soft X-ray (Hudson et al. 1999; Hudson & Schwenn 2000). Temperature 
variation may therefore play a significant role in creating substructures within the cavity. 
However, we note again that the outer boundary of the cavity has similar location and shape 
at all wavelengths (SXR, 171, 195, and 284 A ), indicating that, substructure aside, we are 
observing a true cavity (i.e., a density depletion). 


2.2. Density depletion in cavity 

As noted above, radio and white- light observations have historically been used to demon- 
strate low density in cavities relative to surrounding material. The cavity we are observing 
is rather small compared to those that tend to be studied in white light, which are biassed 
towards sizes easily seen above a coronagraph occulting disk (Gibson et al. 2006). Figure 
5 (a) illustrates this - the white-light streamer is a much bigger structure than the cavity, 
which is visible as a faint depletion near the coronagraph occulting disk (see also Figure 5 
(b)). The location of the cavity within the global white-light structure is indicated by an 
insert showing a density-sensitive line ratio (Fe XII 186/195 A ) obtained by Hinode EIS. 
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Fig. 5. — (a) Large-scale white-light streamer surrounding cavity: MLSO Mk4 white-light (po- 
larized brightness, pB) coronagraph data from August 9, 2007 (19:03:15), and (insert) EIS Fe XII 
186/195 A density-sensitive line ratio from August 9, 2007 (19:46:37). (b) August 9 daily average 
Mk4 r -0 pB map of NE streamer and cavity. The error bar brackets the top of the cavity, with a 
central top height of 1.28 ± 0.0 1 Rq at 43.5 ± 0.5° colatitude, (c) Image showing zoomed-in EIS 
Fe XII 186/195 A line ratio as a greyscale image indicating depleted density in the cavity (red 
contours). EIT 284 A contours and error bar from Mk4 image are overplotted in black. 


This line ratio is related to line-of-sight-averaged density, and is indeed lower in the cavity 
than its surroundings (Figure 5 (c)). Figure 5 (b) shows the cavity in white- light (polarized 
Brightness), which is proportional to electron density integrated along the line of sight, and 
which also indicates depleted density in the cavity. For both of these measurements, some 
sort of assumption about variation of density along the line of sight is required to quantify 
plane-of-sky density. The three-dimensional model that we will establish in this paper will 
serve this purpose (an explicit calculation of plane-of-sky density utilizing it is planned as 
a future analysis, Schmit et al, in preparation). For the moment, we simply make the qual- 
itative note that the depletion of white-light and EIS line-ratio intensity relative to their 
surroundings can only be due to lower density in the cavity, since these measurements have 
little or no temperature dependence. 
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Fig. 6.— 195 A images of cavity from (a) EUVI-B; (b) EIT; (c) EUVI-A. 



Fig. 7. — Relative positions of STEREO-B, Earth, and STEREO-A for 21:00 UT August 9, 2007. 
STEREO-B: east limb Carrington longitude = —37.732°, heliographic latitude = 5.475°. Earth: 
east limb Carrington longitude = —27.575°, heliographic latitude = 6.341°. STEREO-A: east limb 
Carrington longitude = —13.653°, heliographic latitude = 7.111°. 

2.3. Cavity variation with solar longitude 

Figure 1 shows that the cavity was observed for several days, indicating a longitudinally 
extended structure. This can also be seen in Figure 6, which shows observations for August 
9, 2007 from the EUVI instruments on the twin STEREO-A and B spacecraft (21:05:30 UT) 
bracketing SOHO EIT observations (13:12:32 UT) (see Figure 7 for the relative position of 
the STEREO and SOHO spacecraft). These observations confirm that that the structure 
extended a minimum of 24° at that time since the STEREO spacecraft were separated by this 
amount, and the cavity is clearly visible in images from both the A and B viewpoints. Figures 
8 and 9 demonstrate that the actual length is longer, with as much as 90° of longitude showing 
a cavity in Carrington maps, depending on the time and viewing angle of observation. Figure 
1 also indicates that the height of the cavity increases and the center of the cavity at the limb 
moves southward with time, or decreasing longitude. This southward trend is also seen in 
Figures 8 (a) and (c) and in Figure 9; note, however, that it is not apparent in the west limb 
Carrington maps. We will use our forward model to demonstrate how line-of-sight effects 
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Fig. 8.- — Longitude-colatitude Carrington map (CROT2060) using observations of 195 A at the 
solar limb at a height of 1.1 Rq. (a) EUVI B east limb; (b) EUVI B west limb; (c) EUVI A east 
limb; (d) EUVI A west limb. Colored boxes show the location of the cavity observed for each view, 
which differ both in time and viewing angle. Arrows in (a) and (c) show longitude to the left of 
which the EUV cavity contrast diminished. 



Longitude (degrees) 


Fig. 9. — Longitude-colatitude Carrington map (CROT2060) using Nangay Radioheliograph 

432 MHz observations taken from a slice at central meridian at an assumed height of 1.025 R® 
(typical value for this frequency). A data gap occurs near map center. The cavity is seen in the 
orange box; the arrow indicates position to the right of which cavity visibility diminished. 
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Fig. 10. — Aug. 10, 2007 observations of the cavity on the NE limb (indicated by white arrow) by 
Nangay Radioheliograph at 408 MHz (a - b: 11:30:30 UT) and EIT 195 A (c - d: 13:13:46 UT). 
(b) displays radio isocontours on radio image, and (d) displays radio isocontours on EUV image. 
Both data sets have been enhanced for better viewing: radio data with a multi-scale filter, and EIT 
data with a local histogram equalization. The solar limb is shown in (a) and (c) as black line. 

can explain this in Section 4. 

The cavity is the limb manifestation of the filament channel observed on the solar disk 
(see Su et al. (2010) and references therein). This is illustrated by radio observations, which 
show both limb and disk. Figure 10 shows the cavity observations as seen by the Nangay 
Radioheliograph at 408 MHz. The cavity is seen as a brightness depression on the disk 
and extending above the limb to match the extent of the EUV cavity (best seen in the 
isocontours in Figure 10 (b) and (d)). Figure 11 shows the location of magnetic neutral lines 
and filament channels obtained from H-q:. The colored rectangles indicate the Carrington 
longitudes where EUV and radio observations show a cavity, and illustrate that the location 
of the cavity overlies the neutral line of the filament channel observed in H-a, as expected. 
Three of the arrows show the locations where visibility of the cavity is diminished: to the left 
of the yellow and black arrows for the EUV observations, and to the right of the orange arrow 
for the radio observations. The relative visibility of the cavity may vary for these observations 
for several reasons. First of all, the EUV represents line-of-sight integrated limb emission, 
while the radio observations show emission localized to a (frequency-dependent) height at 
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Fig. 11. — H -a synoptic chart (McIntosh 2003) showing magnetic polarity (red-positive, blue- 
negative), neutral lines (black contours), and location of filament channels (green means a filament 
was observed in I I-o) , for Carrington rotation 2060 ±30°. The northern PCFC associated with the 
cavity is clearly visible at all Carrington longitudes. Colored boxes and arrows correspond to those 
in Figures 8 and 9 (note that the red arrow has been replaced with a black arrow for visibility); 
blue arrow indicates top of A-type shape in neutral line (see discussion in Section 4). 

which opacity is near unity. Moreover, the observations represented by the three Carrington 
maps sample three different times: the east limb is observed first by STEREO B and then by 
STEREO A, and then the radio map is based on central meridian slices which are shifted in 
time by several days. Thus time evolution of the cavity may affect its appearance. Finally, 
the Earth- view of the radio observations and the two STEREO spacecraft are from different 
heliolatitudes (or solar B angles), and these differences in combination with the angle and 
bends of the neutral line underlying the cavity can affect its visibility. We will discuss this 
further in Section 4 in the context of the forward modeling. 


3. Three-dimensional cavity morphology 
3.1. Model 

We begin with a three-dimensional, global streamer density model (Gibson et al. 2003). 
This model builds an observationally-constrained corona, beginning with a spherically sym- 
metric background density that matches coronal hole radial density falloff as in Guhathakurta 
et al. (1999). We then superpose a streamer on this background, with a radial density profile 
as defined in Gibson et al. (1999) at the core (centered on longitude/colatitude (<f> 0 , 0 o )). 
Figure 12 shows cartoon representations of the model with streamer parameters labelled, as 
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Fig. 12. — Cartoons illustrating the three-dimensional croissant-like cavity tunnel with elliptical 
cross sections within the Gaussian streamer. Note that for simplicity, the cavity has no tilt so that 
Qctopo = G 0 - (a) view on solar disk, (b) view from above. 

viewed on the solar disk and from above. The streamer density varies in longitude-colatitude 
via Gaussian profiles in its (angular) length and width: the model parameter for the streamer 
base half-length is Si engt h (equivalent to c *^~ 6 ^ in Gibson et al. (2003)), and for the half- width 
is S W idth (equivalent to in Gibson et al. (2003)). The streamer half- width varies lin- 
early from its photospheric value to a “current-sheet” value CS W idth (equivalent to ( - a+b *^ cs ') 
in Gibson et al. (2003)) that remains fixed for all heights greater than i^, simulating the 
finite width streamer stalks surrounding current sheets extending out into the solar wind. 
If the streamer Gaussian has an axis parallel to the solar equator, its length runs parallel 
to longitude and its width to latitude. Otherwise, the length axis (which corresponds to 
an underlying magnetic neutral line) is defined to have an angle m to the equator. The 
streamer may also be non-radial with an angle a to the radial, but, since our cavity is very 
small compared to the streamer (see Figure 5 (a)) any such non-radiality is unlikely to affect 
its appearance. We therefore assume a radial streamer with a — 0. The Gibson et al. (2003) 
model allows multiple streamers, but for our purposes we just need one surrounding the 
cavity that is consistent with white-light observations. Table 1 lists the parameters required 
to define the model streamer, with values corresponding to the fit to be described in the 
next section. 

For this paper, we have modified the streamer model to allow an embedded cavity with 
an elliptical cross-section (see Figure 13 (b)). The cavity is essentially a Gaussian-shaped 
tunnel oriented with its axis along the streamer axis (see Figure 12). The top of the ellipse 
is positioned at (r ct0Po ,9 ct0Po ,(f) 0 ), and the length of the ellipse axis closest to radial ( Cw c ) 


Table 1. Streamer Properties 


Quantity 

Parameter 

Value 

Streamer central colatitude 

0 o 

45.55° ± 0.29 

Streamer central Carrington longitude 

4>o 

-30.33° ± 2.47 

Angle of streamer axis to equator 

m 

-8.95° ± 0.47 

*Tilt of streamer height axis vs. radial 

a 

0° 

*Streamer half-width at photosphere 

S width 

15° 

*Streamer half-length at photosphere 

length 

75° 

* Streamer current sheet height 

Res 

2.5-R© 

*Streamer current sheet half-width 

C S w idth 

1.5° 


Note. — Streamer model parameters; see Gibson et al. (2003) for details. Asterisks indicate 
parameters that have minimal impact on cavity morphology and so are not part of the fit to 
cavity data. Background radial density falloff is defined by the spherically symmetric coronal hole 
background of Guhathakurta et al. (1999) and, for the streamer core, the model of Gibson et al. 
(1999). 


Table 2. Cavity Properties 


Quantity 


Parameter Value 


Cavity top radius at central Carrington longitude (f> 0 

'f'etopo 

1.25 R q ± 0.02 

Cavity top colatitude at (j) 0 

@ctop 0 

47.11° ± 0.70 

Cavity ellipse axis (nearest to radial) at (f> 0 

C'rado 

0.25-R© ± 0.02 

Cavity ellipse axis (normal to C ra d) at (f> 0 

Cnorm 0 

0.19 R q ± 0.01 

Cavity half-length 

C length 

34.5° ± 2. 


Note. — Cavity model parameters. Radial density falloff follows the streamer density falloff, but 
depleted by 50%. 
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Fig. 13. — Elliptical fit to the edge of the cavity as observed in EUV. (a) image from the EUVI-A 
instrument, at 195 A ; (b) we show how an ellipse is fit to user-defined points along the edge of the 
cavity, providing the input needed to determine model parameters. 

is measured down from the top along a line that goes through the streamer core center 
(r = Rq, 0 o , (p 0 ). The bottom of the ellipse is free to lie above or below r — Rq. The other 
ellipse axis has a length C norrrio , defined perpendicular to C ra d 0 at the ellipse center. The 
ellipse size decreases along the streamer length, with a Gaussian half-length Ci eng th required 
to be smaller than Si eng th (see Figure 12 (a)). The bottom of the ellipse remains at a fixed 
radial height along the length of the cavity tunnel, so that the top lowers as the ellipse 
shrinks, maintaining the aspect ratio of the ellipse (Figure 12 (b)). For the purposes of this 
paper, in which we only fit the cavity’s three-dimensional morphology, we assume a radial 
density profile for the cavity scaled to exactly half of the streamer core radial density profile. 

Table 2 lists the model parameters defining the cavity morphology, with values corre- 
sponding to the fit in the next section. The model defined by the parameters then provides 
a fully 3D description of coronal density, which can be integrated along the line of sight to 
reproduce white- light pB data as was done in Gibson et al. (2003). As in that paper, the 
observer’s position, including heliographic colatitude and longitude, can be explicitly taken 
into account in performing line of sight integration. For coronal observables other than white 
light, including EUV and SXR, we must also provide a temperature. For the purposes of 
this paper, we assume an isothermal corona of 1.6 MK for both cavity and streamer. 
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3.2. Determining model parameters from observations 

The model parameters that must be determined, then, are 9 0 ,(f) 0 ,m,r ct0Po ,9 ct0Po ,C ra d o , 
C 'norm 0 -, and Ci engt h (see Tables 1 and 2). In order to fit these, we examine a sequence 
of EUVI images at three-hour time resolution, August 6-12 (EUVI-A) and August 5-11 
(EUVI-B). For each image, we select points along the edge of the cavity, and fit an ellipse to 
these points (Figure 13). As we will describe below, the variation of ellipse vs. Carrington 
longitude (equivalent to time for the rotating Sun) constrains the model parameters. Note 
that although we fit ellipses to both EUVI-A and EUVI-B images, we found that EUVI-B 
images were consistently poorer quality, possibly because STEREO-B was over 13% further 
from the Sun than STEREO-A making the the EUVI-B images dimmer. The parameter 
values quoted in the Tables, and the black lines representing the fits in the figures, thus 
represent a fit to only EUVI-A data. 

The largest source of error is the identification of the ellipse edge, because it is not 
uncommon for nested elliptical structures to appear in a cavity (either internal or external 
to the “true” cavity volume). This may be due to projection effects (e.g., a tunnel with a 
sloping roof might have nested internal ellipses) or to transient brightening of loop- or ellipse- 
like elements internal or external to the ellipse. Generally this ambiguity can be resolved 
by seeking the edge containing an extended (in colatitude) depletion - see, e.g., Figure 3, 
and by examining the cavity over an extended period of time to identify transient effects. 
To quantify the error involved, we have incorporated fits by multiple users (see Figure 14). 
Eight of the co-authors separately performed elliptical fits of the cavity vs. time, resulting 
in a set of model parameters for each person. These parameter values were averaged for 
the final results and are shown along with their standard deviations in Tables 1 and 2. We 
identify the standard error (difference between the multi-user measured values and black line 
model fit, normalized to number of free parameters) within the figure captions. 

The central colatitude of each ellipse-fit at the limb (taken to correspond to the angle of 
the underlying neutral line, 9 n i) is the point where the axis closest to radial ( C ra d ) intersects 
the photosphere. Figure 14 (a) shows 6 n i vs. Carrington longitude, as measured by multiple 
observers, with a linear fit (black line) yielding the angle m of the underlying neutral line 
relative to the solar equator. Figure 14 (b) shows how the maximum altitude of each ellipse 
( fctop ) varies vs. Carrington longitude. The peak in r ctop occurs at Carrington longitude 0 O ; 
given (f) 0 , 9 a follows from Figure 14 (a). 

In our ellipse-fits, the colatitude of the ellipse top {9 (stj0p ) was allowed to vary so that C ra d 
did not have to be truly radial (see Figure 13 (a)). We found that the difference between 
9cto P and 9„i (thus, the degree to which ellipse-fits were non-radial) was small but significant 
(see Figure 15). Table 2 shows the central value 9ct oPo evaluated at <f> 0 for a linear fit of 9 ct0 p 
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Fig. 14. — (a) Cavity central colatitude ( 0 n i ) at r = R@ vs. Carrington longitude, as established by 
ellipse-fits to the cavity as observed in EUVI-A at 195 A . (b) Cavity top height (r ctop ) established 
by same ellipse-fits, also vs. Carrington longitude. Model fit is shown as black line. Different 
color/symbol indicates a different co-author doing the fit (identified by initials at right). The large 
black dot is the location of the cavity top found in white light, see Figure 5 (b). Standard error for 
6 n l is 1.12°, and for r ctop is 0.03 R©. 



Carrington Longitude (deg) 

STEREO-A 

Fig. 15. — Cavity central colatitude 6 n i (black) vs. cavity top colatitude (green). Standard error 
for 6 n i is 1.12°, and for cavity top colatitude is 1.20°. 

vs. Carrington longitude. 

The altitude of the observed cavity bottom was a free parameter in the ellipse-fits. Our 
model assumes a single value for the position of the cavity bottom, which would be best 
observed at longitudes near <f) 0 . For longitudes between (f> 0 ± 15°, we find the cavity bottom 
lay at 1.00.R© ± 0.01. We therefore fixed the cavity bottom at r — Rq for the model, and 
calculated the distance C ra d between the cavity top ( r ctop , 0 rfop ) and the limb (r = Rq, 6 n i) 
for each ellipse-fit. 
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The variation of C ra d vs. Carrington longitude then provides a key constraint, estab- 
lishing a Gaussian variation of the model cavity’s near-radial axis. However, as discussed in 
Gibson et al. (2003), because m 7 ^ 0 the streamer (and so embedded cavity) extends along 
a tilted “heliomagnetic longitude axis”, which is fully determined now by m,9 0 , and 0 O (see 
equations for heliomagnetic longitude (j) mg in section 4.4. of Gibson et al. (2003)). The 
model Gaussian profile for both streamer and cavity are defined relative to this axis rather 
than Carrington longitude. We do a Gaussian fit in this coordinate system (centered on 
0 m£/ — 0 — 0o) which establishes the parameters C ra d 0 (amplitude of Gaussian), and Ci eng th 
(Gaussian half-length along heliomagnetic longitude axis). The central cavity top altitude 
r ctop 0 a l so follows from C ra d 0 and 9 ctoPo , and the assumption that the cavity bottom is at the 
photosphere. Figure 14 (b) shows this fit as a black line plotted vs. Carrington longitude to 
better compare to the data. (Note that we were only able to fit about half of the modeled 
profile of Tctop because of decreasing visibility of the cavity in the EUV images: we discuss 
this further in Section 4.) Figure 16 shows that, although we only used EUVI-A data to 
constrain the fit, it matches EUVI-B observations as well, indicating minimal time evolution 
of the structure. 



Carrington Longitude (deg) 
STEREO-B 



STEREO-B 


Fig. 16. — Points based on EUVI-B data plotted with the EUVI-A model fit (black lines) as shown 
in Figure 14. Standard error for 6 n i is 1.36°, and for rct op is 0.04 We found the quality of data 
was uniformly better for EUVI-A than for EUVI-B, and so used only EUVI-A data in our fit. 


This leaves only one remaining model parameter, C n0 rm o ■ The observed variation of 
ellipse-fit normal axes vs. Carrington longitude somewhat over-constrains the problem, since 
our model requires (heliomagnetic)-longitudinally invariant elliptical cross-sections. This 
essentially forces the same Gaussian profiles for C nor m and C ra d, with only the amplitude of 
the C norm Gaussian to be determined. This is done by fitting 0 mg vs Cnorm with a Gaussian 
and taking the amplitude as C normo . The value of this and the other parameters thus fit are 
recorded in Tables 1 and 2. 
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Fig. 17. — (a) EUVI-A image of the cavity with (b) intensity contours overlaid, (c) Forward 
modeled (line-of-sight integrated) EUV emission using model density and temperature with (d) 
intensity contours overlaid. See online movie for a model-data comparison of the variation of cavity 
location and visibility vs. time/longitude 

4. Results and Discussion 

Figure 17 shows the EUVI-A image of the cavity on August 9, 2007, and the forward 
modeled EUV emission observed from the same Carrington longitude and heliographic lat- 
itude as the STEREO-A spacecraft at that time. The size and shape of the cavity are 
well-matched. For example, we note the enhanced brightness of the poleward rim near the 
limb, which the model reproduces. This asymmetry is a projection effect discussed in Gib- 
son et al. (2003), responsible for the “bananas” observed in white-light Carrington maps 
described in that paper (caused by both foreground and background material projecting to 
higher latitudes). The southward trend of the cavity central colatitude and the variation of 
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Fig. 18. — Forward modeled Carrington map of (line-of-sight integrated) EUV 195 A emission 
from the viewpoints of (a) STEREO A east limb; (b) STEREO A west limb. Maps from the 
STEREO B viewpoints (not plotted) are very similar. Maps show limb emission at the height of 
the cavity midpoint (r = 1.125.R©) 

its visibility with time/longitude can be seen in the online movie associated with this fig- 
ure. This trend translates to an angle of the underlying neutral line relative to the equator, 
illustrated by Figure 18 (a), which shows the forward modeled EUV Carrington maps for 
the STEREO- A east limb viewpoint (compare to Figure 8 (c), bearing in mind we are only 
fitting the NE streamer). 
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Fig. 19. — Carrington map (northern hemisphere) of forward modeled (line-of-sight integrated) 
EUV 195 A emission at cavity midpoint height r = 1 . 1 25 Rq . (a) neutral line angle m and solar B 
angle set to zero; (b) model value of m, but still has zero B angle; (c) model value of m, and the 
solar B angle as seen from the viewpoint of STEREO-A during the time of the cavity observations. 
Note we have shifted to <f> 0 = 180 to center the structure for easier viewing. 

To understand the impact of this angled neutral line on the line-of-sight integrated 
emission, we ran a sequence of forward models. We began with all parameters chosen as in 
Tables 1 and 2, except that the angle of the neutral line relative to the solar equator, m , 
is set to zero and (for simplicity) we set the central longitude (j) Q — 180°. Figure 19 (a) is 
the Carrington map for this model, viewed from a heliomagnetic latitude (solar B-angle) of 
zero. The streamer and its embedded cavity for this case have a mostly symmetric emission 
profile, although the “banana” effect introduces a slight U-shape curvature to the streamer 
and a more obvious enhancement of the northern part of the streamer core (compare upper 


-23- 


vs. lower central red contours). The two vertical red lines bracket the sharply defined part 
of the cavity within the black contour, centered on <j> — <j> 0 . 

In the middle image we run our forward model with the neutral line angle set to the 
model value m — —8.95°, but maintain a heliomagnetic latitude of zero. As discussed in Wiik 
et al. (1994), the optimal viewing angle for a cavity is along an underlying neutral line angled 
at exactly the solar disk polar angle corresponding to its central limb latitude. However, for 
line-of-sight integrated (optically thin) EUV emission, this optimal angle switches sign for 
frontside vs. backside material. In Figure 19 (b), the neutral line angle m tilts the cavity 
axis further away from the solar disk polar angle for the foreground, but more along the 
solar disk polar angle for the background. When the cavity first becomes visible, more of its 
length is in the background, and so the neutral line angle improves the viewing somewhat. 
As we pass the cavity central Carrington longitude 0 O , however, the neutral line angle works 
against visibility. The effect of this is seen in Figure 19 (b), in which the central black 
contour has shrunk and shifted to the right (greater Carrington longitude = earlier times). 

There is one more angle that needs to be taken into consideration: the heliomagnetic 
latitude of the observer. Figure 19 (c) illustrates that, for the STEREO-A viewing position 
(east limb), the solar B angle is oppositely signed and nearly compensates for the neutral line 
angle m. The result is a forward-modeled streamer-cavity structure that is largely symmetric 
in heliomagnetic latitude/longitude, and that is more distinctly visible than the Figure 19 
(b) case. The central black contour is still shifted a little to the right as compared to Figure 
19 (a), but otherwise has a very similar size and shape. 

We now return to Figure 18, which shows our best fit forward modeled EUVI-A emission 
for both for the east and west limb. Figure 18 (a) and (b) can be directly compared to 
Figure 8 (c) and (d), and it provides an explanation for the lack of cavity visibility at the 
west limb. At the west limb, the negative m angle acts in the opposite way from the east 
limb (foreground material tilted down, background material tilted up), and the positive solar 
B angle decreases, rather than increases cavity visibility. Gibson et al. (2003) used a related 
argument to explain the east- west asymmetry in Carrington maps where northern streamers 
slanted up and to the right for east limb maps, and up and to the left for west limb maps 
when viewed from a positive solar B angle (see also discussion in Wang et al. (1997). This 
can in fact be seen in the observations shown in Figure 8, and in the forward model fit of 
Figure 18, where the west limb cavity is not only diminished but is angled up and to the left 
relative to the east limb cavity. 

We were only able to constrain our fit with measurements from ellipses to just past the 
Gaussian midpoint, because, as Figure 8 (c) shows, the cavity contrast was lost below around 
Carrington longitude of 300° (yellow/black arrows). However, despite the loss of contrast, 
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a cavity is observationally distinguishable in Figure 8 (c) down to Carrington longitudes of 
nearly 270° (as indicated by the yellow box). This longitudinal extent is captured by the 
full Gaussian of our forward model (Figure 19 (c)). The loss of cavity contrast may be due 
to departures from a linear m angle of the underlying neutral line. Figure 11 indicates that 
the slope of this neutral line flattens out close to where the EUV cavity dims (yellow/black 
arrows); at this point the slope’s effect on foreground material would be working against 
visibility, and we would lose also the benefit of that slope for the background material. 
Conversely, the A-shaped kink in the neutral line shown by the blue arrow would improve 
visibility, at least for line-of-sight integrated EUV, as there would be less of a harmful slope 
in the foreground, and more of a helpful slope in the foreground. Indeed, the blue arrow 
points to a Carrington longitude where the cavity was particularly visible at the east limb. 


5. Conclusions 

The analysis presented here represents the first step of a process to establish three- 
dimensional distribution of density and temperature in and around a coronal prominence 
cavity. The broad range of observations we have employed, including wavelength and view- 
point coverage, have allowed us to tightly constrain the three-dimensional morphology of this 
cavity. This sort of analysis complements EUV tomographic reconstructions (Vasquez et al. 
2009, 2010), which do not require a predetermined model form for the cavity morphology, 
but which can only provide information for radial heights possessing significant EUV signal. 
Our limb- view ellipse-fitting technique allows us to fit the full three-dimensional cavity mor- 
phology, with the curvature of the cavity boundary providing the clue to where the cavity 
top lies. The comparison with white-light data (e.g., the Mk4 data points in Figures 14 and 
16) allows verification of this extrapolation of shape. 

Perhaps the greatest strength of a “forward-modeling” approach is that it is possible to 
consider how a given three-dimensional structure like the cavity may appear from different 
vantage points and for different observables. The sensitivity of cavity visibility to viewing 
angle is considerable: changes of a few degrees in neutral line orientation and/or observer’s 
solar B angle are enough to obscure it. For example, we have demonstrated with the model 
why a cavity may be clearly visible on one limb and not the other, as the solar B angle can 
act to increase cavity visibility on one limb, but decrease visibility on the opposite limb. It 
is of course possible that evolution of the structure over the two weeks between east and 
west limb viewing may also play a role in the differences. However, the fact that our east 
limb model fit to EUVI-A data reproduced EUVI-B data reasonably well (Figure 16) is 
evidence that the the cavity did not change greatly during the multiple days we observed it. 
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In general, the comparison of forward-modeled observables to the data justifies the form of 
our morphological model. 

Now that we have the three-dimensional morphology of this August 2007 cavity, we 
are positioned to make full use of the unique spectroscopic data from our campaign. The 
data were taken on August 9, which was the day the cavity was most clearly visible from 
Earth. The morphological model presented here allows an explicit inclusion of line-of-sight 
projection effects in interpreting these spectroscopic data: such projection effects may arise 
both because of density and temperature gradients in the cavity and because of the intrusion 
of bright streamer material along the line of sight. Although it is clear that the cavity 
has a relatively low density compared to the surrounding streamer, the quantification of 
its density and temperature remains as a future project of relevance to understanding the 
magnetohydrodynamic state of the prominence cavity. 

Finally, the techniques and software we have developed for fitting cavity observations 
and translating time sequences to a three-dimensional morphological model represents a 
considerable asset for future analyses of cavities. The techniques are general, and may be 
applied to cavities observed at a variety of wavelengths. It will be interesting to use them 
to study how cavity properties may vary with wavelength, size, latitude, and time of solar 
cycle, as well as how they may evolve just prior to an eruption. 
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